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riFORMANCE OF CONVOLUTION CODING CONCATENATED 


WITH MFSK MODULATION IN A GAUSSIAN CHANNEL 


1. Introduction 

In a m-ary frequencj^- shift-keyed communication link (MFSK) one of the M 
messages, ; (f = 1, 2, 3, ... M) is transmitted during each T second interval 
as a sinusoidal tone at frequency f£ . This tone is then modulated onto a carrier 
frequency f^ to produce the single transmitted frequency f^ + f^ . Ideally, at 
the receiver the carrier frequency is removed and the resulting signal fre- 
quency f£ is detected using a spectral analysis receiver. In order to reduce the 
probability of a bit error and to increase the available bit rate for the same 
system parameters (error rate, transmitter power, range etc.) concatenation of 
convolution coding with coded MFSK communication link has been proposed. In 
this report we shall study the improvement in db due to concatenation over the 
conventional m-ary coding. The performance of convolution coded system de- 
pends on the channel parameters R [ 1 ] and E„(p) [21. For a discrete 
memoi’y less channel with M inputs and J outputs. 


E (p) = max 
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j = 1 k= 1 

where {pj) are the a priori input probabilitic ;, and 
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It can be shown by random coding arguments that the undetected error 
probability using Fano decoding algorithm satisfies the following bounds: 


r A2~'' ""comp/V for <R 

J N comp 

I k2-'^P\ for R > R 

’ N comp 


where A is a constant in the order of unity. 


K is the code constraint length. R^j is the code rate in bits per transmitted 
waveform, p is the pareto exponent. The pareto exponent p is the solution of 
the equation 


Eq(P) 

p 




( 3 ) 


The decoder computational distribution must be known to determine the 
necessary buffer size. If c is a random variable equal to the number of de- 
coder computations required to decode an information bit, then c has a pareto 
distribution, i.e., 

Prob (c > L) = BL-^ (4) 

where B is a constant in the order of unity. It follows from eq. (4) that p must 
be greater than 1 for finite average computation. Therefore Eo(l) = Rcomp 
called the computational cut off rate. 


A space communication channel can be a .^urately modelled on an additive 
Gaussian noise channel. Digital data is transmitted over this channel by coding 
it into a set of analog waveforms suitable for transmission over the channel. A 
demodulator converts the received signal back into digital form. From the point 


2 



of view of the encoder and decoder the channel consists of the combination of 
the modulator, white Gaussian noise channel, and demodulator. The relative 
efficiency of various modulation and demodulation schemes can be compared by 
calculating Ej,/Nq necessary for R^omp “ where is the signal energy 
per information bit and N q is the one-sided noise spectral density. In this 
report the results of calculations for orthogonal modulation with non-coherent 
detection and Q-level correlator quantization are presented. 

2. Channel Model 

It will be assumed that the signal transmitted is one of the M equal energy, 
orthogonal signals 


Sj (f ) = e. (t) sin j 

= 1 , 2 , . . . M 

( 5 ) 


s,(t) s.(t)dt =E^S.,, 

“'0 

83. = 1, Y3 = i 

= 0, ^ i 

(6) 


The received signal is 

y(t) = (t) sin (copt + 6>) + n(t) (7) 

where 0 is a random variable uniformly distributed over (0, 2 t 7), n(t) is white 
Gaussian noise with spectral density N^/2. 

The maximum likelihood receiver for this non-coherent channel calculates 
the M quantities 

Yi = l/xf + ojf , for i = 1 , 2 , . . . M (8) 


3 



where 




rT 


y(t) e. (t) sin ojo(t) dt 


( 9 ) 


'■ =lP^ 

y N E 


rr 


y(t) ej (t) cos c^o(t) dt 


( 10 ) 


The device that performs these calculations is known as an envelope detector. 

The correlator outputs are statistically independent and have probability densities 


^2 2 , 
(y; + a ) 


f(yi/Si) = yi e 


io(®yi) 'J(yi) 


= fs.n(yi)> j " i 
2 


(11) 


f(yi/Sj ) = yj e 2 u(yp 


= f„(yj )> for j i 


( 12 ) 


where 


a 



u( ■ ) is the unit step function. 


In this study the correlator outputs are quantized to one of the Q levels. The 
receiver output is a vector consisting of a list of the M correlator quantum 
levels. Therefore the channel has Q“ possible outputs and M possible inputs. 
With this method the optimum output can be approached by increasing fine 
quantization. 
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3. R for Unquantized Outputs 

comp *■ 

Rcomp for exact optimum correlator output case can be found by letting 
the quantization become infinitely fine. From eqs. (1) and (2) it follows that 


R 


c omp 


max 





(13) 


where Y is the correlator output vector and f(Y| s\) is the joint density for the 
correlator outputs given input Sj. For symmetric channels is a maxi- 

mum for equally likely inputs, i.e., Pj = 1/M for all i . On substituting eqs. (11)- 
(12) in eq. (13) it follows that 


M 


^ CO 

2 

1 + (M - 1) 

(ax)dx 
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f^comp vs Ef^/Ng for M = 2, 4, 8, 16, 32, 64 is plotted in the enclosed figures. 

4. Rco^p with Q-Level Quantization 

To approximate the exact output case the half line (0, o3) can be divided into 
a disjoint set of intervals Aj such that A . = (0, ®). If the correlator output 
y^ ^ j , then the quantized value is taken to be yg = j. Let 


P(yi=i/se)= I t.,„(x)dx = p^^„(i) 
■'Ai 


for f = j 


= I f„(x),x =p„(i: 

i 


for f / j 
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From eqs. (1) and (2) it follows that 


R = log, 

comp *=* 4 


M 


1 + (M-1) 


Z /Ps.n(i)Pn(i) 


i = 1 


(14) 


For simplicity equal width quantization intervals were chosen. The quan- 
tizer input-output relationship is illustrated below. Optimum thresholds were 
searched for by trial and error. The optimum thresholds depends on the signal- 
to-noise ratio. Therefore, sub optimum thresholds good for a wide range of 
signal-to-noise ratios (0-19) were chosen and shown in Table 1. Rcomp vs E^/Nq 
for these thresholds are plotted. 
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5. Eq(p) VS p with Q-level Quantization. 

For symmetric channels equally likely mputs, i.e., p. = l/M for all i , 
maximize the expression for Eq(p). UsL:g the fact that the correlator outputs 
are statistically independent, equation (1) reduces to 
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n i+p 

using the sub optimum thresholds of Table 1 calculated values of Eo(/0) vs P for 
various values of En/Nq and M are plotted. 


M M 

z: L ^ 
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